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Abstract 



It is known that neutrino oscillations may map <p e : 4>^ : 4> T = 1 : 2 : 0, 
the initial flavor ratios of ultrahigh-energy neutrino fluxes produced from a 
distant astrophysical source, into : : (ff = 1 : 1 : 1 at the detector of 
a neutrino telescope. We remark that this naive expectation is only valid in 
the fjL-T symmetry limit, in which two neutrino mixing angles satisfy # 13 = 
and # 2 3 = 7r/4. Allowing for the slight breaking of \i-r symmetry, we find 
0° : 0° : 0° = (1 - 2A) : (1 + A) : (1 + A) with A characterizing the 
combined effect of # 13 ^ and 23 ^ ir/4. Current neutrino oscillation data 
indicate —0.1 < A < +0.1. We also look at the possibility to probe A by 
detecting the V e flux of E- ~ 6.3 PeV via the Glashow resonance channel 
V e e — > W~ — > anything. 
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[X] The solar [1], atmospheric [2], reactor [3] and accelerator [4] neutrino experiments 
have convinced us of the existence of neutrino oscillations and opened a new window to 
physics beyond the standard model. Given the basis in which the flavor eigenstates of 
charged leptons are identified with their mass eigenstates, the phenomenon of neutrino 
mixing can simply be described by a 3 x 3 unitary matrix V which links the neutrino flavor 
eigenstates (v e , v , v T ) to the neutrino mass eigenstates (i^, u 2 , u 3 ): 

(1) 

In the so-called standard parametrization of V [5], V e2 = sin# 12 cos9 13 , V e3 = sm9 13 e~ tS 
and V^ 3 = sin 9 23 cos 9 13 . Here we have omitted the Majorana CP-violating phases from V, 
because they are irrelevant to the properties of neutrino oscillations to be discussed. A global 
analysis of current experimental data (see, e.g., Ref. [6]) points to 9 13 = and 9 23 = n/4, 
which motivate a number of authors to consider the \x-r permutation symmetry for model 
building [7]. In the limit of \x-r symmetry, the (effective) neutrino mass matrix M v takes 
the form 



M v = V M U V T = lb c d , (2) 
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where M v = Diag{m 1? m 2 , m 3 } with m i (for i = 1,2,3) being three neutrino masses, and 
V is a special pattern of V with 9 13 = and 9 23 = tt/A. Note that 9 12 is arbitrary and 
5 is not well-defined in V . The parameters a, b, c and d of M u are in general complex 
and can produce the desired Majorana CP-violating phases for V . Any slight breaking of 
fi-r symmetry in M v will result in non- vanishing 9 13 and small departure of 9 23 from 7r/4. 
Nontrivial 5 can also be generated from the breaking of fi-r symmetry, leading to the effect 
of CP violation in neutrino oscillations. 

The main purpose of this paper is to investigate how the effect of fi-r symmetry breaking 
can show up at a neutrino telescope. We anticipate that IceCube [8] and other second- 
generation neutrino telescopes [9] are able to detect the fluxes of ultrahigh-energy v e (]7 e ), 
v (V ) and v T (u T ) neutrinos generated from very distant astrophysical sources. For most of 
the currently-envisaged sources of ultrahigh-energy neutrinos [10], a general and canonical 
expectation is that the initial neutrino fluxes are produced via the decay of pions created 
from pp or pj collisions and their flavor content can be expressed as 

K> <t>r} = {\, f, O}0 O , (3) 

where <f) a (for a — e, r) denotes the sum of v a and V a fluxes, and O = <p e + M + <p T 
is the total flux of neutrinos and antineutrinos of all flavors. Due to neutrino oscillations, 
the flavor composition of such cosmic neutrino fluxes to be measured at the detector of a 
neutrino telescope has been expected to be [11] 

K, * 0?} = {\, \, \}<i>o- (4) 
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However, it is worth remarking that this naive expectation is only true in the limit of /i-r 
symmetry (or equivalently, 9 l3 = and 9 23 = vr/4). Starting from the hypothesis given in 
Eq. (3) and allowing for the slight breaking of \x-r symmetry, we are going to show that 

{€■■ = {(1-2A):(1 + A):(1 + A)} (5) 

holds to an excellent degree of accuracy, where A characterizes the effect of /i-t symmetry 
breaking (i.e., the combined effect of 9 13 ^ and 9 23 ^ 7r/4). We obtain —0.1 < A < +0.1 
from current neutrino oscillation data. We find that it is also possible to probe A by 
detecting the V e flux of E v m 6.3 PeV via the well-known Glashow resonance (GR) channel 
z/ e e — > W~ — > anything [12] at a neutrino telescope. 



Let us define 0( D ) = 0f, D ) + 0- D ' ) (for a — e, a, r) throughout this paper, where 0( D ) 

and (f)^ denote the v a and V a fluxes, respectively. As for the ultrahigh-energy neutrino fluxes 
produced from the pion-muon decay chain with <f> = 0-^ = 0, the relationship between <f> 
(or <f) v ) and <#? (or <$ ) is given by D = P ea + P or 05 = P ea + P 



m 



which Pg Q and Pg a (for ct = e, /i, r and /3 = e or //) stand respectively for the oscillation 
probabilities -P(^ — > f a ) and i^X 7 ^ — > X 7 ^). Because the Galactic distances far exceed 
the observed neutrino oscillation lengths, Pp a and Pp a are actually averaged over many 



oscillations. Then we obtain Pg Q = P^ a and 



^ = EI^H^I 2 , (6) 



i=l 



where and (for a, (3 — e, /i,t and i — 1,2, 3) denote the matrix elements of V defined 
in Eq. (1). The relationship between a and 0° turns out to be 

€ = <PeP ea + 4>,P,a ■ (7) 

To be explicit, we have 

0e = f (Pee + 2Pj , 

# = | + 2P MT ) ■ (8) 

It is then possible to evaluate the relative magnitudes of 0^, 0^ and 0^ by using Eqs. (1), 
(6) and (8). 

In order to clearly show the effect of /i-r symmetry breaking on the neutrino fluxes to 
be detected at neutrino telescopes, we define a small quantity 

b = 9 23 -^, (9) 

Namely, e measures the slight departure of 6*23 from 7r/4. Using small 9 13 and e, we may 
express \V ai \ 2 (for a = e, //, r and i — 1,2, 3) as follows: 
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(10) 



Combining Eqs. (6) and (10) allows us to calculate Pg Q . After a straightforward calculation, 
we obtain 

P ee + 2P^ e = l-e sin 2 29 12 + ^ sin4£ 12 cos 5 + 0{e 2 ) + 0(9 2 13 ) , 
P efl + 2P^ = l + £ - sin 2 29 12 - °f sin4£ 12 cos 5 + 0(e 2 ) + 0(£ 2 3 ) , 



P er + 2P^ T = l + £ - sin 2 2£ 12 - °f sin 4£ 12 cos 5 + 0(e 2 ) + 0(9 2 13 ) . 
Substituting Eq. (11) into Eq. (8), we arrive at 



(11) 



f (1-2A) , 



# = ?(1 + A) , 



f (1 + A) , 



where 



A = i (2£sin 2 2^ 12 - 13 sin 40 12 cos + C(e 2 ) + 0(9 2 13 ) 



(12) 
(13) 

holds. Some 



Eq. (5) is therefore proved by Eq. (12). One can see that 
discussions are in order. 

(1) The small parameter A characterizes the overall effect of fi-r symmetry breaking. 
Allowing 5 to vary between and n, we may easily obtain the lower and upper bounds of 
A for given values of 9 12 (< tt/4), 9 13 and e: -A bound < A < +A bound , where 



Abound = \ (2|e| sin 2 29 12 + 9 13 S inA9 12 ) + 0{e 2 ) + 0(9 2 3 ) 



(14) 



It is obvious that A = — A bound when e < and 5 = 0, and A = +A bound when e > and 
5 = it. A global analysis of current neutrino oscillation data [6] indicates 30° < 9 12 < 38°, 
9 13 < 10° (« 0.17) and |e| < 9° (« 0.16) at the 99% confidence level, but the CP-violating 
phase S is entirely unrestricted. Using these constraints, we analyze the allowed range of A 
and its dependence on 5. The maximal value of A bound (i.e., A bound ps 0.098) appears when 
\e\ and 9 13 approach their respective upper limits and 9 12 ~ 33° holds, as one can clearly 
see from Fig. 1(A). Indeed, we find that A bound is not very sensitive to the variation of 9 12 
in its allowed region. 
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Provided 9 13 = holds, we easily obtain A bound = 0.5|e| sin 2 29 12 < 0.074 when 9 12 
approaches its upper limit. If e = (i.e., 6*23 = 7r/4) holds, nevertheless, we find A bound = 
0.25# 13 sin4# 12 < 0.038 as 9 12 approaches its lower limit. We observe that A bound is more 
sensitive to the deviation of 9 23 from tc/4. 

(2) Of course, A = exactly holds when 9 l3 = e = is taken. Because the sign of e and 
the range of 5 are both unknown, we are now unable to rule out the nontrivial possibility 
A rs in the presence of 9 13 7^ and e 7^ 0. In other words, A may be vanishing or 
extremely small if its two leading terms cancel each other. It is straightforward to arrive at 
A pa from Eq. (13), if the condition 

— = cot2# 12 cos5 (15) 

is satisfied. Due to | cos5| < 1, Eq. (15) imposes a strong constraint on the magnitude of 
s/9 l3 . The dependence of e/9 l3 on 5 is illustrated in Fig. 1(B), where 9 12 varies in its allowed 
range. One can see that \e\/9 13 < 0.6 is necessary to hold, such that a large cancellation 
between two leading terms of A is possible to take place. It should be remarked again that 
the above result is a natural consequence of the assumption made in Eq. (3) for the initial 
flavor ratios of ultrahigh-energy neutrino fluxes. 

The implication of Fig. 1 on high-energy neutrino telescopes is two-fold. On the one 
hand, an observable signal of A 7^ at a neutrino telescope implies the existence of significant 
/i-r symmetry breaking. If a signal of A 7^ does not show up at a neutrino telescope, on 
the other hand, one cannot conclude that the \x-r symmetry is an exact or almost exact 
symmetry. It is therefore meaningful to consider the complementarity between neutrino 
telescopes and terrestrial neutrino oscillation experiments [13], in order to finally pin down 
the parameters of neutrino mixing and leptonic CP violation. 

(3) To illustrate, we define three neutrino flux ratios 



R e = 



0° 



R 



+ ' 



* s «TSj?' (16) 

which may serve as the working observables at neutrino telescopes [14]. At least, R can be 
extracted from the ratio of muon tracks to showers at IceCube [8] , even if those electron and 
tau events cannot be disentangled. Taking account of Eq. (12), we approximately obtain 

R e « - - -A , 

e 2 2 ' 

1 3 . 

Rr « l + fA . (17) 

It turns out that R e is most sensitive to the effect of \i-r symmetry breaking. 
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As a straightforward consequence of = 0° shown in Eq. (12), R = R T holds no 
matter whether A vanishes or not. This interesting observation implies that the "/i-r" 
symmetry between R and R T is actually insensitive to the breaking of \i-t symmetry in the 
neutrino mass matrix M v . If both R e and R are measured at a neutrino telescope, one can 
then extract the information about A from the difference of these two observables: 

R,-R e = \^- (18) 

Taking A = A bound » 0.1, we get R^-R e < 0.22. 



We proceed to discuss the possibility to probe the breaking of \i-r symmetry by 
detecting the V e flux from distant astrophysical sources through the Glashow resonance (GR) 
channel V e e — > W~ — > anything [12]. The latter can take place over a narrow energy interval 
around the V e energy E^ R « M^ / /2m e ~ 6.3 PeV. A neutrino telescope may measure both 
the GR-mediated 17 e events (A^ R ) and the v^+V ^ events of charged-current (CC) interactions 
(iV"CC _ ) j n the vicinity of E9 R . Their ratio, defined as R RG = N^ R /N^ C +V , can be related 
to the ratio of u^s to u^s and u^s entering the detector, 



J5 

v. 



R " - Ktw t , ' (19) 

Note that 0? , 0° and 05^ stand respectively for the fluxes of 77 e 's, u^s and 77^'s before the 
RG and CC interactions occur at the detector. In Ref. [15], -Rqr = °-^o w hh a 30.5 has 
been obtained by considering the muon events with contained vertices [16] in a water- or 
ice-based detector. An accurate calculation of a is certainly crucial for a specific neutrino 
telescope to detect the GR reaction rate, but it is beyond the scope of this work. Instead, 
here we concentrate on the possible effect of \i-r symmetry breaking on R . 

Provided the initial neutrino fluxes are produced via the decay of 7r + 's and 7r~'s created 
from high-energy pp collisions, their flavor composition can be expressed in a more detailed 
way as follows: 

K, K K> *»*> M = {I > \> h I ' ' °} 00 • (20) 

In comparison, the flavor content of ultrahigh-energy neutrino fluxes produced from p~f 
collisions reads 1 

K, *v <^ V tf„ T , ^} = {\, , I , 1 , , O}0 O . (21) 



Note that the dominant reaction to generate electron and muon neutrinos in p"f collisions is 
j>7 — > A + — > 7r + n with ir + — > [i^v^ and p + — > e + v e V IJl . There is no production of V e , because the 
produced neutrons can escape the source before decaying [17]. In contrast, the numbers of i/ e 's 
and F e 's in Eq. (20) are identical as a result of the equal amount of 7r + 's and 7r~'s produced from 
inelastic pp collisions. 
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For either Eq. (20) or Eq. (21), the sum of 0„ and <j) v (for a = e, //, r) is consistent with 
<f) a in Eq. (3). 

Due to neutrino oscillations, the V e flux at the detector of a neutrino telescope is given 
by 0° = 0- P ee + 0- P^ e . With the help of Eqs. (6), (10), (20) and (21), we explicitly 
obtain 

$SPP) = y(l"2A) , 

^.(P7) = ^(sin 2 2^ 12 -4A) . (22) 

The sum of 0° and 0° ; which is defined as 0° has been given in Eq. (12). It is then 
straightforward to calculate R by using Eq. (19) for two different astrophysical sources: 

1 3 
R (PP) ~ 2 ~ 2 A ' 

„, , sin 2 26» 12 4 + sin 2 26» 12A 

^o(P7) ~ — ^ i ^ • (23) 

This result indicates that the dependence of R (pp) on 9 12 is hidden in A and suppressed by 
the smallness of 9 l3 and e. In addition, the deviation of R Q (pp) from 1/2 can be as large as 
1.5A bound pa 0.15. It is obvious that the ratio R Q (jry) is very sensitive to the value of sin 2 29 12 . 
A measurement of R (py) at IceCube and other second-generation neutrino telescopes may 
therefore probe the solar neutrino mixing angle 9 12 [15]. Indeed, the dominant production 
mechanism for ultrahigh-energy neutrinos at Active Galactic Nuclei (AGNs) and Gamma 
Ray Bursts (GRBs) is expected to be the p~y process in a tenuous or radiation-dominated 
environment [18]. If this expectation is true, the observation of R (py) may also provide us 
with useful information on the breaking of fi-r symmetry. 



In summary, we have discussed why and how the second-generation neutrino tele- 
scopes can serve as a striking probe of broken ji-T symmetry. Based on the conventional 
mechanism for ultrahigh-energy neutrino production at a distant astrophysical source and 
the standard picture of neutrino oscillations, we have shown that the flavor composition 
of cosmic neutrino fluxes at a terrestrial detector may deviate from the naive expectation 
0D ; 0D . 0D = 1 . l . L Instead; 0D : 0D : 0D = ^ _ 2 A) : (1 + A) : (1 + A) holds, where 

A characterizes the effect of /j,-t symmetry breaking. The latter is actually a reflection of 
9 13 7^ and 9 23 ^ 7r/4 in the 3x3 neutrino mixing matrix. We have examined the sensitivity 
of A to the deviation of 9 13 from zero and to the departure of 9 23 from 7r/4, and obtained 
—0.1 < A < +0.1 from current neutrino oscillation data. We find that it is also possible 
to probe the breaking of ji-T symmetry by detecting the V e flux of E- pa 6.3 PeV via the 
Glashow resonance channel V e e — > W~ — > anything. 

This work, different from the previous ones (see, e.g., Refs. [13-15,19]) in studying how to 
determine or constrain one or two of three neutrino mixing angles and the Dirac CP- violating 
phase with neutrino telescopes, reveals the combined effect of 9 13 ^ , 9 23 ^ 7r/4 and 5 ^ n/2 
which can show up at the detector. Even if A ^ is established from the measurement of 
ultrahigh-energy neutrino fluxes, the understanding of this \x-r symmetry breaking signal 
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requires more precise information about 8 l3 , 6 23 and 8. Hence it makes sense to look at 
the complementary roles played by neutrino telescopes and terrestrial neutrino oscillation 
experiments (e.g., the reactor experiments to pin down the smallest neutrino mixing angle 
9 n and the neutrino factories or superbeam facilities to measure the CP-violating phase 5) 
in the era of precision measurements. 

The feasibility of our idea depends on the assumption that we have correctly understood 
the production mechanism of cosmic neutrinos from a distant astrophysical source (i.e., via 
pp and p7 collisions) with little uncertainties. It is also dependent upon the assumption that 
the error bars associated with the measurement of relevant neutrino fluxes or their ratios 
are much smaller than A. The latter is certainly a challenge to the sensitivity or precision 
of IceCube and other neutrino telescopes under construction or under consideration, unless 
the effect of fi-r symmetry breaking is unexpectedly large. Nevertheless, there is no doubt 
that any constraint on A to be obtained from neutrino telescopes will be greatly useful 
in diagnosing the astrophysical sources and in understanding the properties of neutrinos 
themselves. Much more efforts are therefore needed to make in this direction. 

The author would like to thank J.X. Lu for warm hospitality at the Interdisciplinary 
Center for Theoretical Study of USTC, where part of this paper was written. He is also 
grateful to H.B. Yu and S. Zhou for useful discussions. This work is supported by the 
National Natural Science Foundation of China. 
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FIG. 1. (A) the dependence of A bound on 6 12 , where 9 13 and e take their respective upper limits 
(i.e., 6 13 < 10° and e < 9°); (B) the nontrivial condition for A = 0, where 6 12 is allowed to vary in 
the range 30° < 12 < 38°. 
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